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NOTATION 



English 

a speed of sound, fect/second 

skin friction coefficient, dimensionless 

c specific heat at constant pressure, 

1 foot -pounds/slug- degree Rankine 

F dimensionless ratio of various boundary 

layer thicknesses 

h specific enthalpy, foot-pounds/slug 

k shear layer mixing coefficient, dimensionless 

L body length, feet 

M Mach number, dimensionless 

m mass flux, slugs/second-square foot 

n normal coordinate, feet; also 

velocity profile integer 

p static pressure, pounds/square- foot 

Q heat added per unit mass, foot-pounds/slug; 

also dimensionless heat addition rate of Eqn. (12) 

q heat addition rate, foot-pounds/slug -second 

R universal gas constant, foot-pounds/slug-degree 

Rankine 

Re Reynolds number, dimensionless 

r radial coordinate, feet 

S specific entropy, foot -pounds /slug -degree Rankine 
s streamline coordinate, feet 

T static temperature, degrees Rankine 

U velocity vector, feet /second 

U,u speed, feet/ second 
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w speed along streamline, fect/second; also 

nondimens ion ali zed by stagnation speed 

of sound, a 
’ s 

w/ r dimensionless heat zone length 

z/r dimensionless position on body or shear layer, 
measured from base 

Greek 

a flow field type, dimensionless 

3 (M 2 - l) h 

y ratio of specific heats, dimensionless 

6 boundary layer thickness, feet 

<5* boundary layer displacement thickness, feet 

6** boundary layer momentum thickness, feet 

e dimensionless small quantity, 0 < e << 1 

n left running characteristic coordinate, feet; 

also dimensionless heating parameter of Ref. 7 

G flow angle, radians 

k dimensionless ratio of average internal to 
external speed 

y Mach angle, radians 

v Prandtl -Meyer angle, radians 

£ right running characteristic coordinate, feet 
p density, slugs/cubic foot 

<f> perturbation velocity potential, square feet/second 
?2 curl U, radians/second 

00 angular velocity, radians/second 

Superscripts 

( )' dimensionless perturbation quantity 
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Subscripts 

1 average condition 

00 freestream condition 
b base condition 

e condition at edge of shear layer 

L laminar 

o,s total or stagnation condition 
T turbulent 

tr laminar transition point 
trt turbulent transition point 
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I. INTRODUCTION 



Projectiles in supersonic flight experience drag. One 
source of drag is due to the pressure acting on the base 
(rear surface) of the projectile, which, under normal cir- 
cumstances, is below ambient pressure. The base drag force 
is equal to the base surface area normal to the direction 
of flight multiplied by the mean difference between the 
ambient and base pressures. 

External burning has been studied by many authors [Ref. 

1] to determine its ability to reduce the drag of a body or 
even to produce thrust. The basic objective of external 
burning in a supersonic flow, Fig. 1, is to create compres- 
sion waves which impinge on the base flow region and increase 
the static pressure within the recirculation zone, thereby 
increasing the pressure acting on the base of the body to 
reduce base drag or produce thrust. 

Other methods of reducing base drag or producing thrust 
are available. The rocket assisted projectile provides thrust 
from combustion of self-contained fuel and oxidizer. With 
external burning, additional fuel can replace some of this 
oxidizer since oxygen is used from the atmosphere for burning 
the excess fuel. The ram jet assisted projectile^ delivers 
excellent performance but does require air inlets and lias 
additional hardware costs. The method oi base burning offers 
advantages similar to external burning, but has not proven 

its ability to create thrust. 

12 
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Figure 1. External Burning Projectile 



The external burning assisted projectile provides a 501 
increase in range when fuel rich propellant is substituted 
for one third of the conventional projectile payload. In 
addition to an increase in range, external burning can be 
used to decrease the projectile flight time to the target, 
increase its kinetic energy at the target and provide a 
flatter trajectory. 

The base pressure ratio, p^/p TC , is a measure of the ef- 
fectiveness of the external burning process. As shown in 
Fig. 2, the behavior and effect of this ratio can be related 
to the flow angle, 0, of the boundary layer streamline at 
the base. 

For purposes of analysis, the external burning process 
may be described thermodynamically as a heat addition term 
in the energy equation, bearing in mind that numerous other 
factors must be considered with actual combustion, such as 
the effect of pressure on combustion rate. The linearized 
problem involving two dimensional planar uniform flow with 
unit heat addition has been solved by Tsienand Beilock in 
closed form [Ref. 2,3]. In supersonic flow it was shown 
that a pressure perturbation, due to the heat source, is 
transmitted along a characteristic. Further study of the 
linearized flow equations by Rues has shown that there is an 
equivalence between heat, force and mass sources [Ref. 4). 
These unit sources were given an integral formulation by 
Fuhs to determine the flow field perturbation resulting from 
a source zone [Ref. 5] . It is the study of these linearized 
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results which suggested that for experimental investigation 
the external burning process may be simulated by supersonic 
nozzle wall contours which generate identical characteristics. 

Heat addition in a two dimensional planar supersonic flow 
was investigated by Strahlc, using one dimensional and lin- 
earized heat addition theory, and shown to reduce base drag 
or produce .thrust [Ref. 6,7]. In this study, the base flow 
was described by the Crocco-Lees [Ref. 8] model, and the 
resulting base pressure ratio was determined for heat addi- 
tion zones of different intensity, location and extent. Ex- 
perimental investigation [Ref. 9-11] has demonstrated the 
feasibility of external burning and shown that it can in- 
crease base pressure. Heat addition in an axisymmctric 
supersonic flow has been studied by Hosack and O'Leary using 
the method of characteristics and the Korst base flow model 
[Ref. 11] . 

The object of the present investigation was to study 
base pressure in two dimensional planar and axisymmetric 
supersonic flow with external heat addition using the method 
of characteristics and the Crocco-Lees base flow model. The 
effects of varying Reynolds number, including boundary layer 
transition, and projectile spin have been included. A con- 
current experimental investigation studied the effect on 
base pressure of various axisymmetric nozzle contours which 
simulated the external burning process. 

It should be mentioned that several combinations of ex- 
ternal heat addition analysis and base flow models have not 
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been investigated theoretically. These include two dimen- 
sional planar method of characteristics and linearized theory 
with the Chapman-Kors t base flow model, and linearized axi- 
symmetric flow with both the Chapman-Korst and Crocco-Lees 
base flow models. The linearized axisymmetric flow problem 
will be discussed later. 
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II. THEORY 



The theory employed for this investigation must describe 
the steady flow in two entirely different flow regimes, 

Fig. 1. The external flow .is supersonic and can be accurately 
treated by the method of characteristics. The internal vis- 
cous flow, or recirculation zone, is dissipative and diffi- 
cult to describe. As in the two dimensional case, it was 
assumed that the Crocco-Lees model [Ref. 1,12-15] adequately 
describes the axisymmetric internal flow. Boundary layer 
type (laminar or turbulent) , thickness and transition were 
included as a function of Reynolds numer since these effects 
profoundly influence the base flow behavior. 

The characteristic equations have been derived for flow 
with heat addition [Ref. 16,17]. In vector notation the 
steady equations for flow without viscosity, heat conduction 



or body forces arc 




V- (pU) = 0 


(1) 


pU-ViJ + Vp = 0 


(2) 


U • V (h + %U-U) = q 


(3) 


p = pRT 


(4) 



where p is density, U is velocity, p is static pressure, h is 
enthalpy, q is the time rate of heat addition per unit mass, 

R is the universal gas constant and T is static temperature. 
Using a natural coordinate system [Ref. 18], Fig. 3, the 
characteristic equations were derived from the above equations 
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Systems 
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in the following convenient form (Appendix A) 
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where 0, v and y are respectively the flow, Prandtl-Meyer 
and Mach angles, U is speed, S is specific entropy and h^ 
is total enthalpy, h^ = h + %U 2 . The left and right running 
characteristics are n and £ respectively. The natural co- 
ordinates, s and n, are along and normal to the streamline. 

For two dimensional planar flow, a is 0, and a is 1 for 
axisymmetric flow. The characteristic equations, (5) and 
(6) , have been multiplied by the characteristic length of 
base radius, r^ , to obtain a dimensionless form. Therefore, 
in contrast to the dimensional coordinates n, £ and n used 
in Fig. 3, the coordinates n , £ and n in Eqns. (5) and (6) 
are understood to be dimensionless. Also, the specific heat 
at constant pressure, c , was assumed constant in the enthalpy 

relation h = c T. 

P 

The characteristic computations with Eqns. (5) and (6) 

1 9 S 

were simplified by neglecting the term {j(T-^ - — ^ )/i' b , 
which is the vorticity. One may assume the external flow to 
be irrotational before heat addition and then determine the 
restrictions on heat addition which preserve irrotati onali ty . 
Crocco's theorem for steady flow may be written (Appendix B) 



TVS + UXfi = Vh 



0 



( 7 ) 



where 0 = VXU. Taking the curl of Eqn . (7) 
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U-Vft = ft-VU - ftV*U + VTXVS . 



( 8 ) 



Assume ft = 0 before and during heat addition. Then, from 
Eqn . (8), the condition which preserves irrotationality is 

VTXVS =0. (9) 

If one redefines the heat addition and has uniform upstream 
conditions at infinity [Ref. 19], then 

V(h Q - Q) = 0 (10) 

where Q is the heat per unit mass which has been supplied to 
the fluid element. Using Eqn. (7), with IT = 0 , Eqns . (9) 

and (10) yield the following result 

VQXV(hU 2 ) = 0 (11) 

which is the condition preserving irrotationality; that is, 
the gradients of kinetic energy and heat added must be par- 
allel. If Q is a small perturbation, (Q/h OT ) = e, then Eqn. 
(11) is of order e 2 and the irrotationality condition is 
automatically satisfied to order e. To neglect the vorticity 
term in Eqns. (5) and (6), one must insure 



r b q _ U • VQ 

Uc T r bUc T 
P P 



e << 1. 



( 12 ) 



In the computations discussed in Section V, the quantity 
of Eqn. (12), (r^q) / (Uc T) , was assigned a constant value 
independent of local pressure, within the heat zone and was 
taken as zero outside the heat zone. Thus, within the heat 
zone, VQ was approximately parallel to the gradient of kinetic 
energy and the irrotationality condition of Eqn. (11) was 
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satisfied. This condition was not satisfied along the inner 
and outer streamlines of the heat zone, since a jump in Q 
existed normal to the streamlines. However, by assigning a 
value for Eqn. (12), (r^q)/(UCpT) < 0.2 << 1, the irrota- 
tionality condition was satisfied to order epsilon along the 
streamline boundaries of the heat zone. 

Equations (1) through (4) may be combined to give 



V-U kj*V( yj*U) = -AfT 

a 2 v J cl 

P 



(13) 



where a is the speed of sound. Assuming small perturbation, 
Eqn. (13) may be linearized (Appendix G) 

V 2 4> - M* 0 = q' (14) 

where M m is the Mach number at infinity, q' is the heat ad- 
dition, q = h^q', <(> is the perturbation velocity potential, 

U = iU ro + U' , U ' = Vcf) and T = T to ( 1+T'). In Cartesian co- 
ordinates, Eqn. (14) may be written 

- 3 2 <}> + (j> + (p = q' (15) 

T xx ^yy T zz H K - 

where 3 2 = - 1. For supersonic flow, the solution of 

Eqn. (15) is found by integrating the sources in the fore- 
cone volume, [Ref. 20], 



4> = - 



_JL 

2 7T 



1 

V r 



q ' dxi dy i dz i 

((x-x,) 2 - 6 2 ((y-y i ) 2 * (z-z,) 2 TP5 



(16) 



Although the linearized theory was not pursued further, it 
should be noted that this theory could be used instead of 
the method of characteristics to establish the effect of 
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external burning at the base flow shear layer. Studies of 
linearized supersonic flow through ducts have provided solu- 
tions for an annular source distribution [Ref. 21-23], which 
could be used in conjunction with centerline sources and 
sinks [Ref. 18] to describe the axisymmetric external flow 
with burning. 

The two dimensional planar Crocco-Lees theory [Ref. 8] 
considers the base flow as a separated boundary layer merg- 
ing into a free jet profile and derives integral equations 
to describe the mean values of the flow quantities (Appendix 
C) . From the integrals of the conservation equations, a 
single ordinary differential equation is derived. This or- 
dinary differential equation has a saddle point singularity; 
existence of the singularity is a key feature of the Crocco- 
Lees model. The behavior of this equation is partly deter- 
mined by an empirical F - k relation and by the nature of 
the shear layer mixing. k is the ratio of average internal 
base flow velocity to the external velocity, u^/u^. F is a 
ratio formed by using the various boundary layer thicknesses. 
The turbulent shear layer mixing coefficient, k r} . , was as- 
signed the constant value, 0.03, by Crocco-Lees. The laminar 
value used for this investigation was 0.0033, about 1/9 the 
turbulent value. 

F is conceptually a profile form factor; when the veloc- 
ity profile changes shape, the value of F changes. As dis- 
cussed in reference 8, the F - k functional relationship is 
based on many known viscous flows. For example the Blasius 



solution is represented by a point in the F - k plane. 
Falkner-Skan boundary layer solutions, which involve either 
favorable or adverse pressure gradients, become a curve in 
the F - k plane. The velocity profiles of jets under cer- 
tain conditions are described by the error function; these 
jets become a point in the F - k plane. Likewise the Howarth 
boundary layer solutions for flows in pressure gradients are 
a curve in the F - k plane. 

The F - k curve corresponds to viscous flow near walls, 
i.e. boundary layers, as well as free jets. Furthermore, 
adverse, zero or favorable pressure gradients are represented 
by positions on the F - k curve. Heat addition in a super- 
sonic flow generates compression waves. When z/r < 0, the 
compression waves impinge ahead of the base; the boundary 
layer on the body and the shear layer experience an adverse 
pressure gradient. When z/r > 0, the compression waves im- 
pinge downstream of the base and only the shear layer flows 
into the adverse pressure gradient established by the exter- 
nal burning. Since the F - k curve incorporates pressure 
gradients, the presence of an adverse pressure gradient due 
to heat addition should not alter its validity. As stated 
in reference 8, the F - k curve needs additional experimental 
verification; this applies to the case with external burning. 
Further discussion of the F - k curve may be found in refer- 
ences 6, 12, 13, 14 and 15. 

The reformulation of the Crocco-Lees theory for axisym- 
metric flow [Ref. 1,12-14] resulted in the same single ordinary 



24 



differential equation which had been derived for planar flow. 
In two dimensional planar flow this equation is numerically 
integrated by starting at the singularity, or so called 
critical point, and working toward the base, using the 
Prandtl -Meyer relation to obtain flow angle for various Mach 
numbers. With external heat addition, one may no longer use 
the Prandtl -Meyer relation for flow angle at the shear layer, 
since the heat addition, described by either linearized or 
one dimensional flow, introduces a correction term [Ref. 6,7]. 

This procedure cannot be used for the axisymmctric prob- 
lem without heat addition due to the 1/r term arising in the 
continuity equation; thus, one must start at the base and 
integrate toward the critical point. As was anticipated, 
attempts to numerically integrate through this saddle point 
singularity were unsuccessful; however, by repeated integra- 
tion one may determine two Mach numbers on the downstream 
side of the expansion fam which straddle the saddle point 
singularity, thus establishing the base pressure within the 
limits given by the two Mach numbers, Fig. 4. This technique 
was used with Fqns . (5) and (6) for both axisymmctric and 

two dimensional planar flow. 

To incorporate the effects of various boundary layers on 
the body (laminar, transitional, turbulent) and shear layer 
transition, it was necessary to formulate a transition model, 
Fig. 5. While this model treats the boundary layer mass flux 
differently, it is not a departure from the existing Crocco- 
Lces theory. In fact, this transition model extends the 

Crocco-Lces base flow model. 
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Figure 5. Transition Model 



■Since the transition from laminar to turbulent flow is 
more or less probabilistic in nature and, depending on the 
Reynolds number, may occur on the body or along the shear 
layer, it was assumed that the transition process could be 
described by an integrated normal distribution and two 
transition Reynolds numbers. The transition Reynolds num- 
bers .identified the start and end of the transition process, 
and the integrated normal distribution ranged from 0.02 to 
0.98 in the span between the two specified transition 
Reynolds numbers. This assumption replaces the previously 
assumed transition curves [Ref. 8]. 

The transition model must provide several parameters 
which are required in the Crocco-Lees base flow formulation. 
It must establish the values for the mixing coefficient, k, 
along the shear layer and the boundary layer mass flux at 
the base to permit determination of F. 

The boundary layer was assumed to be at constant static 
pressure and total enthalpy, and the velocity profile was 
assumed to have a 1/n dependence, where n is an integer, 

1 < n < 7. This formulation may be integrated (Appendix D) 
for both planar and axisymmetric flows to give the boundary 
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completely turbulent, the n is taken as 7 and 6 is computed 
from the compressible one-seventh power solution [Ref 8], 

6 t = (0.037L/(fcc b ) 1/5 ) (C f /C £ ) (6/6**) . The term (C fM /C f .) 
is a ratio of skin friction coefficients for compressible 
and incompressible turbulent flow [Ref. 24] . It was recog- 
nized that the boundary layer thickness computed for planar 
flow will be different from the axisymmetric thickness. 

This difference will be small if 6/r^ and radial velocity, 
v/U, are of order e. The dimensional term pv/r from the 
axisymmetric continuity equation will be of order c ? upon 
nondimensionali zing with <5 and U. Thus, the governing equa- 
tions and their solutions are identical to order e in both 
cases . 

The first transition Reynolds number, Re £r , is where the 

laminar flow starts transition to turbulent. The second 

transition Reynolds number, Re £rt , is where the flow must be 

completely turbulent (Re £ £ = 4Re £r ) . As one changes values 

of base Reynolds number from Re £r to ^ e tr£ > the values for n 

and 6, which determine the boundary layer mass flux at the 

base, must vary from 1 to 7 and to 6 £ , respectively. This 

transition was computed by the error function, erf, since it 

is readily available in FORTRAN and can be directly related 

to the integral of the normal distribution. A numerical 

span, S, of 4, which takes the value of erf from -0.96 to 

0.96, was used to normalize the difference (Re. . - Re. ). 

v trt tr' 



Let Re. = J*(Re. + Re„ . ) 
tra v tr trt' 

S 



Then 



ERF B - CRe tra • R3 b )/(R0 trt - Rc tr )). (17) 
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So‘, when Re^ > Re tr 



6 



BASE 



= s l ch + erf b ) + 6 t (% - ERF fi ) 



(18) 



and 



n BASE 



= Hk + erf b ) + 7 c% - erf b ) 



(19) 



were the values used to compute the boundary layer mass flux 
at the base (Appendix F) . 

When transition from laminar to turbulent flow occurred 
on the shear layer, the value for the mixing coefficient, 
k, at a particular point was computed in a similar fashion 
using the error function to transition from laminar mixing, 
kj , to turbulent mixing, krj, , Fig. 6. 

The effect of shear layer mixing on base pressure is 
shown schematically in the momentum control volume of Fig. 7. 
The axial pressure on surfaces cd and de , and momentum flux 
leaving the volume, de , must be balanced by the pressure on 
surfaces ab and be and momentum flux entering the volume, 
abed. Thus, to maintain this balance, an increase of momen- 
tum flux along cd due to mixing requires a corresponding de- 
crease of base pressure, p^ , along ab . 

One additional assumption was necessary to obtain good 
agreement with experimental data. At low Reynolds numbers 
and completely laminar flow, it was necessary to increase 
the value of the laminar mixing coefficient, kj , in direct 
proportion to the increase of the thickness of the boundary 
layer at the base. This assumption may be expressed as 



k L * k ( Re tr /Re b ),S 



( 20 ) 
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Figure 6. Transition of Shear Layer Mixing Coefficient 
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Figure 7. Momentum Control Volume Depicting Shear Layer Mixin 



where Re tr is the transition Reynolds number, Re^ is the 

Reynolds number of the body, based on length or chord, L, 

and k is the laminar mixing coefficient when Re 4 = Re, . 

tr b 

Some justification may be provided for Eqn. (20) by consider- 
ing the Blasius boundary layer solution where the mass flux, 

_ P 

m, is given by 3. 28pux/ (Re) 2 and the Crocco-Lees theory 

where kpu is taken equal to dm/dx. If one forms the deri- 

P — t 

vative of 3 , 28pux/ (Re) 2 , dm/dx is found to vary as l/(Re) 2 , 

thereby suggesting the form assumed in Eqn.' (20). Although 
k^, was held constant at 0.03 for all calculations, a similar 
analysis using the incompressible turbulent solution suggests 
that k^. might vary as l/(Re)'*‘^'\ 

It has been suggested that the effect of projectile spin 
could be included in the mixing coefficient [Ref. 15], k, 
by using k = k(l+rw/U), where rw/U is the spin number. 

However, since k . always increased with an increase in 
spin number, and, since it has been determined through numer- 
ical experiments with the Crocco-Lees theory, Fig. 8, that 
increasing k decreases the base pressure, this simple hypoth- 
esis contradicts experiment [Ref. 25,26] where, in some 
cases, increasing the spin number increased the base pressure. 
In Fig. 9 it is clear that at certain values of Reynolds num- 
ber, Re^ , an increase of Reynolds number resulted in an in- 
crease in base pressure. This fact suggests that the spin 
number could be used to shift the p^/p^ versus Re^ curve, 

Fig. 9, instead of increasing the mixing coefficient. Sev- 
eral methods were attempted at this point, such as increasing 
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gure 9. Computed Values of p, / p vs. Reynolds Number, Mach 2.0, No Heat Addition 



Re^ 5y [Ref. 27] (l+(^-jj) 2 )~ or (l+(^j) 2 ), but the method 

which provided the best correlation was to simply decrease 
the transition Reynolds numbers. The spin transition 
Reynolds numbers used to compute p^/p^ versus Re^ were form- 
ed by dividing the corresponding no spin transition Reynolds 
numbers by (l+(^j)). Interpretation of base pressure data 
for a spinning projectile is given in Section V. 
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III. EXPERIMENTAL INVESTIGATION 



To carry out experiments related to the five inch pro- 
jectile, the mid point of a maximum range trajectory was 
selected. This flight condition established an altitude 
of 23,000 feet and a Mach number of 2.0 for determining 
similarity parameters. For the five inch projectile, this 
gave a body Reynolds number of 15.7x10^, a radius Reynolds 
number of 1.51x10^, and a turbulent boundary layer thick- 
ness to radius, 6/r, of 0.1313. The spin at 23,000 feet 
was assumed to be the same as the muzzle spin, 250 revolu- 
tions per second, which gave a spin number, rw/U, of 0.1595. 
Using these nomial values, it was found that all similarity 
parameters could be matched with a free jet Mach 2.0 coaxial 
nozzle exhausting at room pressure if the center cylinder 
had a radius of 0.641 inches, 6 of 0.084 inches, and a spin 
of 800 revolutions per second. As stated in the introduc- 
tion, the actual combustion process was simulated with 
various nozzle contours which duplicated the compression 
waves generated by the external burning. 

The simulation process is depicted in Fig. 10 where the 
streamline deflected by the presence of a heat zone can be 
replaced by a nozzle wall to generate an identical flow 
field . 

The free jet experimental facility, Figs. 11 and 12, was 
conventional in most respects and has been previously re- 
ported fRef. 28,29]. The use of cast aluminum filled epoxy 
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Nozzle Wall 



Replaces Streamline 




Flow Field With 
Compression 



(b) Identical Flow Field With Compression Generated by 
Nozzle Wall to Simulate Heat Zone 



Figure 10* Heat Zone Simulation With Nozzle Wall 
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'Figure 12. Schematic Diagram of Free Jet 
Experimental Facility 
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norzlcs, Fig. 13, was the only new hardware innovation. The 
technique of casting employed a wax plug contoured to the 
desired shape with an aluminum template. This plug was 
turned slightly undersize and then restored to size using a 
thin coat of automotive body plastic. After melting the wax 
from the cast nozzle, the plastic peeled away from the epoxy 
to leave a smooth nozzle wall. Without this plastic, the 
heat generated by the epoxy had softened the wax and caused 
some irregularity in the nozzle wall. 

The location of the zone of combustion to be simulated 
by the nozzle was estimated using an empirical formula for 
penetration deptli of a jet normal to a supersonic frec- 
stream [Ref. 30], For Mach 2.0 flow, internal pressures 
from 300 to 500 psia were computed to provide penetration 
to 2% inches with various orifice diameters. Thus it was 
assumed that external burning at 5 inches from the projec- 
tile centerline would be practical. Assuming a fuel rich 
propellant exhaust, ARC-168, with an available heating value 
of 7200 BTU/lb and 0.6 lb/sec flow rate, it was found using 
one dimensional theory that the flow should be deflected 
about 4h° for a distance of 15 inches, or three diameters, 
downstream. This computation used the integral methods of 
reference 5 which gave a perturbation velocity normal to the 
flow of -163 feet per second. The resulting flow angle 
tangent was -163/2058, which gave a computed flow angle of 
-4.54 degrees along the inner streamline of the heat zone. 
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Figure 13. Coaxial Epoxy 



Nozzle . 



Insert Showing Spacers. 
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These estimated values for the heat zone were taken as 
the upper limit and were used to form a test matrix and for 
subsequent nozzle design. Several different nozzle contours 
were computed by the method of characteristics for 4?i° stream- 
line deflection at a radius equal to a cylinder diameter 
from the centerline and with different assumed lengths for 
the combustion zone. The test matrix included using spacers, 
Fig. 13 insert, with each nozzle so that the first full com- 
pression waves impinged ahead of, at, or behind the base of 
the cylinder. Pressures were measured along the base and 
side of a stationary centerline cylinder, Fig. 13, using 
mercury manometers. Attempts to carry out experiments with 
a spinning cylinder were unsuccessful due to faulty air bear- 
ing design (Appendix I). Time did not permit correction of 
this difficulty. 

The differences in the flow field with the stationary 
cylinder and different nozzles can be seen in Fig. 14. The 
photograph of the uniform flow nozzle in Fig. 14a used a 
horizontal knife edge and the expansion fan at the base ap- 
pears dark. The base pressure ratio was 0.635. In Fig. 14b 
the vertical knife edge was positioned to make a positive 
density gradient from the centerline to the right appear 
brighter. The flow with compression, Fig. 14c, had an ex- 
pansion fan at the base but achieved a higher base pressure 
ratio of 1.065 due to compression upstream of the base as 
well as along the shear layer. The compression upstream of 
the base increased the pressure right before the expansion 
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Figure 14a. 
Uniform Flow, 
Horizontal 
Knife Edge 



Figure 14b. 
Uniform Flow, 
Vertical 
Knife Edge 



Figure 14c. 
Compression Flow, 
Horizontal 
Knife Edge 



Figure 14. Schlieren Photographs at Mach 2.0 Showing 
Uniform Flow and Flow With Compression 
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fan. so that it was greater than p^, thus the flow was near 
ambient pressure after ■ passing through the expansion fan. 

Some test results using the stationary cylinder are given 
in Table I. Except for the nozzle with w/r = 4.34, these 
tests indicate that compression applied upstream of the 
base, z/r = -0.5, is less efficient at increasing base pres- 
sure than compression applied at the base, or slightly down- 
stream. The nozzle with w/r = 4.34 had a strong shock wave 
at the exit, due to compression waves reflected from the 
cylinder [Ref. 15], which, with z/r = -0.5, impinged on the 
base flow region ahead of the critical point. In turn, in- 
creased pressure across the shock wave was impressed on the 
subsonic flow in the recirculating zone and gave a higher 
base pressure. The shock wave observed with this experimen- 
tal nozzle would not be expected with uniform external 
burning . 

In Table I the experimental base pressure ratios, p^/p^ , 
are shown as the upper number in each box for the various 
test positions, z/r, of each nozzle, w/r. The cylinder 
radius is r, the axial location where the first full strength 
compression wave impinges on the body or shear layer is z 
(z=0 at the base) and w is the length of the heat zone. The 
corresponding Reynolds number, Re^ , is shown below each value 
of z/r. The lower number shown in each box is the computed 
specific impulse for the assumed projectile flight condition 
at 23,000 feet. The uniform flow nozzle, with w/r = 0.0, 
was taken as the reference pressure ratio to compute the 
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Table I. Base Pressure Ratio and Specific Impulse 
With Simulated Heat Addition 
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base pressure force created by the other nozzles. The fuel 
flow rate was taken proportional to (w/r)/(6.0), since 6.0 
was the combustion zone length assumed for the fuel flow 
of 0.6 lb/sec. 

Once Mach number and operating pressure and temperature 
are selected, the remaining variable controlling Reynolds 
number is the length or scale of the apparatus. The free 
jet facility of Figure 11 has sufficient mass flow to per- 
mit duplication of Reynolds number at 23,000 feet altitude 
provided the cylinder has sufficient length. Consequently 
the nominal Reynolds number was chosen as 15.7X10^ for 
23,000 feet to match capability of the facility. 

When the design of the spinning cylinder was initiated 
it became apparent that a long cylinder had problems in 
regard to shaft whip or the first vibrational mode in 
bending. A shorter cylinder reduced the Reynolds numbers 
to the values shown in Table I. Conditions at a higher 
altitude were simulated, or else conditions for a short 
projectile at 23,000 feet were duplicated. 

Two cylinders could have been designed. Each cylinder 
would have required a separate family of nozzles. A choice 
was made to use a short cylinder for both stationary and 
spinning tests. 

Variation in z/r as shown in the sketch with Table I 
changes Reynolds number slightly. The variation in Re^ in 
Table I is due to this fact. 
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IV. NUMERICAL INVESTIGATION 



The computer program for this numerical investigation was 
written to calculate either two-dimensional planar or axisym- 
metric problems. The program used the method of characteris- 
tics where the axisymmetric term in Eqns . (5) and (6), 

r^ (sinpsinO) /r , was deleted for two-dimensional planar flow. 
Similarly, the effect of radius in the mass flux calculation, 
at the base and along the shear layer, was deleted or re- 
tained as required. The characteristic mesh, Fig. 15, was 
computed from the body or shear layer upstream through the 
heat zone, if present, to a reference freestream Mach line 
(i.e., backwards along a right running £ characteristic). 

The error from the reference condition was used to correct v 
on the body, since 0 was given, or 0 along the shear layer 
since v was computed from the Crocco-Lees theory. Each body 
or shear layer point was iterated until the error at the up- 
stream boundary condition was less than 10 radians. At 
each point in the characteristic mesh, the value of p was 
determined for the resulting value of v by the method given 
iii Appendix H. In retrospect, this method of characteristics 
computation could be made more efficient if one could devise 
a scheme which works with arbitrary points along the body or 
shear layer and increments the characteristic mesh along the 
reference freestream Mach line. Then only one calculation 
would be required along each characteristic, rather that the 
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Figure 15. Characteristic Mesh 



iterative scheme used in the present method. This approach 
was considered, but it was not clear how one would formulate 
the calculation for the corner or the shear layer mass in- 
flux (Appendix E) . 

The flow field was calculated along the body to the base 
at which point all computed values along the last character- 
istic were stored. These values were used to reset the flow 
field for each iteration of the straddle method. 

The Crocco-Lees theory was formulated in terms of w £ 
instead of Mach number, where w is the velocity at the edge 
of the shear layer divided by the freestream stagnation speed 
of sound, w = U/a . To start the straddle method a value 
of w was chosen for the downstream side of the expansion 
fan which was too large. The expansion fan was divided into 
ten parts; downstream of the fan the computation proceeded 
toward the critical point. Each new point along the shear 
layer, computed by the Crocco-Lees theory (Appendixes C, E 
and E) , required several iterations along the characteristic 
to the upstream Mach line boundary condition, where the flow 
angle error after each iteration was used to correct the 
shear layer flow angle. 

Integration in the F - w g plane followed the first path 
shown in Fig. 16. It can be seen that this path falls to 
the right of the path passing through the saddle point 
singularity, x, and, as the value of F was decreased, the 
slope, dF/dw , became negative at 1'. This occurred when 
the denominator used to form dF/dw became positive [Ref. S] , 
i.e. both the numerator and denominator of dF/dw were 
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negative from 1 to 1'. This transition from positive to 
negative slope with positive denominator occurred when the 
initial guess for w c was too large. By systematically de- 
creasing the initial guess for w g , say by tenths, one 
eventually has the integration path from 2 to 2', which is 
similar to the previous case, except that the negative 
slope at 2' is associated with the numerator, used to form 
dF/dw c , becoming positive. Thus, the last two initial 
guesses for w g have straddled the saddle point singularity. 

The interval between these two values may be divided, say 
in hundredths, and the process repeated to find two values 
of w c which more closely straddle the saddle point singularity, 

The heat zone was specified to start and stop at certain 
values of axial position, z, and to lie between an inner and 
outer streamline. The quantity r^q/Uc^T was assigned a con- 
stant numerical value, between 0.15 and 0.175, within the 
zone, which gave the desired -4^° deflection to the inner 
streamline. External to the heat zone r^q/Uc T was zero. 

Before deflection, the inner and outer streamlines have 
values of approximately 2.0 and 3.0, maximum body radius, 
respectively. These values were chosen to match the rate of 
compression increase used in the nozzle design. As the com- 
putation proceeded downstream, the values of radius and flow 
angle for the inner and outer streamlines were computed to 
insure realistic boundaries for the heat zone. 

When a sufficient amount of properly located heat was 
added to the external flow, the base pressure ratio was 
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computed to be greater than 1. It is interesting to note 
that for this condition the method of characteristics rou- 
tine generated a compression fan at the base, rather than 
the usual expansion fan. The compression fan calculation 
mathematically folds the flow field on itself [Ref. 31] . 

This of course is physically impossible and is replaced in 
real flow by an oblique shock. For values of p^/p^ = 1.5 
at Mach 2.0, the loss of total pressure due to the oblique 
shock is less than II, which corresponds to an entropy in- 
crease, AS/R, of 0.01. This resulted in an error of less 
than 0.18° in the turning angle, which was neglected. There 
is also the possibility of boundary layer separation up- 
stream of the base for values of p^/p^ greater than one 
[Ref. 7] . This effect was not considered in the calcula- 
tions. The computer program is discussed in Appendix J. 
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V. DISCUSSION OF RESULTS 



The computed values of P^/p OT versus Re^ , shown in Fig. 

9, were based on the apparent transition Reynolds numbers 

from the uniform flow nozzle experiments. These were 

Re^ r = 1.95X10^ and Re ^ = 7.8X10^. The curves display 

the same general features and approximate values found 

from experiment [Ref. 32-35]. It should be noted in Fig. 9 

that for both the two dimensional planar and axisymmetric 

flow the maximum value of p^/p^ occurred below Re and 

the minimum value of p,/p occurred below Re. . . 

1 b 1 00 trt 

A slightly different choice of transition Reynolds 

numbers. Re. = 1.5X10^ and Re. . = 6.0X10^ resulted in 
7 tr trt 7 

the curve of p^/p,^ versus Reynolds number shown in Fig. 17. 

One may compare the curves of p^/p^ in Fig. 9 and Fig. 17 
with the experimental curves shown in Fig. 18. While the 
curves are by no means identical, they do agree in having 
the same general shape. Figure 19 shows additional experi- 
mental data to compare with Fig. 9. The locations of maximum and 
minimum computed p^/p^ are controlled by the choice of values 

for Re. and Re. . . The ratio of h for Re. /Re. . used in 
tr trt tr trt 

these computations was based on information in reference 8. 

This ratio is an approximation which varies considerably for 
different bodies, depending on their shape and surface con- 
dition. The maximum value of p^/p^ was controlled by the 
choice of kj , while the choice of k.^ controlled the second 
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Reynolds Number 

Figure 18. Experimental Values of P b /p vs. Reynolds Number Without Heat Addition, 
trom Ref D 32. w , 
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Figure 19. Experimental Values of p T /p vs. Reynolds Number Without Heat Addition 



maximum value of p^/p^ to the right of Re t « The slight de- 
crease of P^/Poo to the right of this second maximum is in 
agreement with the r/6 sensitivity curve of Fig. 8, since 
increasing Re decreases 6 and therefore increases r/6. 

Three additional curves are shown in Fig. 17; these are 
transition and critical point locations. The laminar and 
turbulent transition points are found simply by dividing 
base Reynolds number into Re tr and Re^ t . The critical 
point location was taken from calculation with the Crocco- 
Lees equations along the shear layer where dF/dw g changed 
signs; it is the singularity in Fig. 16 which was discussed 
in Section IV. In going from small to large values of 
Reynolds number, the shear layer between the critical point 
and the base gradually changes from completely laminar to 
completely turbulent as the laminar and turbulent transition 
points move upstream to the base and then onto the body. 

Where this occurs, the shear layer mixing coefficient, k, 
gradually changes from k^ to k^,. Increasing from some small 
value of Reynolds number, the thickness of the laminar 
boundary layer at the base decreases until the laminar 
transition point reaches the base and moves onto the body. 

As the boundary layer at the base transitions to turbulent 
its thickness increases. The thickness continues to increase 
as the turbulent transition point moves forward along the 
body, reaching a maximum when the transition point is near 
the nose of the body. Beyond this point, the thickness 
decreases . 



4 

In Fig. 17, the value of P^/Peo starts to decrease slight- 
ly before reaching Rc^ , while the laminar transition point 
is still on the shear layer and the base boundary layer is 
laminar. Although the boundary layer thickness is still 
decreasing at this point and in turn should provide an in- 
crease in p, /p^, the value of k along part of the shear 
layer (between the laminar transition point and the critical 
point) starts to increase from k^ to k,p and consequently 
dominates the overall change in P^/p TO with Reynolds number. 

The mechanism responsible for the abrupt minimum value 
of p^/p^ in Fig. 18 was not clear. However, this minimum 
was calculated and is shown in Fig. 17 as a dashed line. 

From the sensitivity curve in Fig. 8, it appeared that k 
would be most likely quantity responsible. In numerical 
experiments this abrupt minimum was calculated by assuming 
additional mixing. The additional mixing was assumed to 

have a normal distribution about i.e. k = k^,(l+ 

“ x ^ / 2 ^ 

(e ' )/(2tt) 2 ) . The argument, x, was adjusted so that the 

additional mixing fell to about ten percent of its maximum 

value in 2X10^ either side of Re. . . This additional mixing 

trt 

was not used in any other calculations since it was not 
understood and was not clearly observed in the experimental 
data reported later in this section. 

Calculations were made for external heat addition in 
two dimensional planar flow to compare with Strahlc's 
previously reported calculations [Ref. 7]. The values of 
Rc ti , = 1.95X10^ and Re^^ = 7.8X10^ were assumed and a 
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Reynolds number of Re^ = 20.0X10° gave the same value of 
6/r used by Strahle. The Mach number was M = 2.3. The heat 
addition zone was started at the same point, z/r = 0.4, and 
assigned the same length, w/r = 3.0, but to avoid subsonic 
flow with the method of characteristics it was necessary to 
increase the height of the heat zone from the value of 0.2r^ 
used by Strahle to a value of one base radius. The values 
of Q = r^q/UCpT were reduced accordingly to provide the 
same total heat addition based on the rectangular heat zone 
used by Strahle, i.e. Q = 0.2ri/tany. The results in Table 
II show good agreement at the low intensity of heat addi- 
tion, but only fair agreement at the higher intensities with 
the method of characteristics program showing higher values 
of P^/Poo- Although the non linearity included in the method 
of characteristics calculation may be responsible for part 
of the higher base pressure ratio, the main cause was at- 
tributed to expansion of the heat zone. This effect was 
included in the bounding streamline calculation for the heat 
zone in the method of characteristics program and resulted 
in a greater amount of total heat added than in the linear- 
ized theory calculations. 
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Table II. Comparison of Calculated Base Pressure Ratio with 
External Heat Addition, Two Dimensional Planar Plow and Mach 
Number of 2.3. 
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Both experimental and calculated results are shown in 
Fig. 20 for a Mach number of 2.0. The value of p^/p^, = 
0.6314 at Rc^ = 10.4 was the average of nine runs with the 
uniform flow nozzle. The probable error about this average 
was ± 0.004, with a maximum observed deviation of - 0.0104. 
In comparing the experimental values of Fig. 20 and Table I 
with early test results using the free jet, variations as 
high as 7% in P^/Poo were found for similar runs. These 
variations were attributed to problems, such as mechanical 
alignment during nozzle installation or leaks in the mano- 
meter tubing, which were found and corrected during subse- 
quent testing. 

Early test data were rejected and were not reported 
herein. As experience was gained with operation of the 
facility, the 7 % variations were narrowed to less than . 

Cylinder wall temperatures were also recorded during 
these nine runs using a thermocouple imbedded in the wall, 
1H inches from the base'. These temperatures ranged from 
5°F to 35°F. A plot of P^/p ro versus temperature showed 
no correlation scatter in p^/p^ due to temperature was less 
than 1.41. Since the stagnation temperatures of the air 
supply were not recorded, they could not be compared with 
the cylinder wall temperatures. 

The values of Re. = 1.95X10^ and Re. . = 7.8X10^ were 

t r t r t 

assumed for the computed values of p^/p^ shown in Fig. 20 
as triangles. The experimental values of p^/p ro for the 
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uniform flow nozzle without heat addition showed good agree- 
ment with the computed values. Three of the four experimen- 
tal values of p^/p^, with w/r = 3.91 compression agreed with 
the computed values. The observed value of p^/p^ with 
z/r = -0.5 was less than the computed value. This differ- 
ence was attributed to an adverse pressure gradient along 
the cylinder wall that was not considered in the calcula- 
tion. In view of this result, base pressure ratios computed 
by the program given in Appendix J should be considered with 
caution if the full compression due to external burning im- 
pinges ahead of the base on the projectile body. 

Figure 21 shows a comparison of experimental [Ref. 26] 
and computed values of base pressure ratio resulting from 
spin . 

In the Crocco-Lces model there are several input param- 
eters including M to , k^, , Re , ^ e trt anc * ^ " K relation- 

ship. A procedure for predicting the influence of spin on 
P 5 /P 00 was discussed previously at the end of Section II. 

This procedure was used for the computed values in Fig. 21. 

Choice of the input parameters remains. Preliminary 
calculations were conducted with M =2.86, Re. = 1.95X10^ 
and = 4.6X10^. Values of p^/p^ had the correct trends 

relative to experimental values. However, the calculated 
P b /Pco ^ or ^ e b = 2.6X10^ and zero spin was about 0.54. To 
decrease the calculated P^/P TO > the Mach number was arbitrar- 
ily increased to 3.0. For Re^ = 2.6X10^, zero spin and M to 
= 3.0, the calculated and observed p^/p TO values agree. In 

6 3 



0.5* -t- 



0.45 



Pb 

P 



oo 



0 X 



20 X 



30 X 



40 X 



♦0 



*20 

*30 

«40 



2 6 

X Experimental Value' 1 , 

Mach 2.86 

* Computed Value, Mach 3.0 
Re tr = 1.95X10 6 , Re trt = 4.6X10 6 

Spin: 0, 20,000 RPM, 30,00ORPM 
and 40,000 RPM 



.0 



0.4 .. 




20 X 



#30 

#40 



30 X 



.40 

'0 

*20,30 



HU « 

30 . 
20 . 



0 * 



0.35- 



2.6 



40 x 



0,40 X 
20,30 X 



40 X 

30 X 
0,20 X 



3.0 3.4 

Reynolds Number, million 



3.8 



2 6 

Figure 21. Comparison of Experimental and Computed 
Values of Base Pressure Ratio vs. Reynolds 
Number With Spin 



64 



retrbspect, additional calculations should have been con- 
ducted with increased mixing, i.e. larger k.^., while retain- 
ing correct = 2.86. In fact, spin should increase k.j, 
as discussed in Section IT. 

A better procedure for calculating the influence of spin 
might be to use a shift in transition Reynolds number coup- 
led with a k = k (1+ (rm/U) ) , as discussed by Smithey, 

Naber, Caswell and Fuhs [Ref. 15] . The dashed segment of 
Pp/Poo vs • curve i n Pig- 17 shows the influence of in- 

creased mixing (larger k^,) in the vicinity of minimum p^/p^. 

The computed values shown in Pig. 21 required four hours 
time on the IBM 360/67. In view of the computing time, the 
procedure suggested in the preceding paragraph was not 
accomplished . 

To provide some estimate of the optimum location for 
external heat addition, a fixed fuel flow was assumed to be 
burned in zones of different length and location. The ini- 
tial radius of the inner streamline was 2.0 and the outer 
streamline initial radius was 3.0. The heat zone were 
assigned integer values from 1 through 6 for w/r. The zone 
with w/r = 1.0 was assigned a value of Q = r.q/Uc T = 0.2. 

D p 

This value was decreased by l/(w/r) for the other five zones. 
The heat zones were positioned along the axial axis starting 
with compression at the base and moving downstream in in- 
teger values of radius. The results of two of these compu- 
tations are shown as the two lower curves in Pig. 22. The 
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Figure 22. Computed Values of Base Pressure Patio at 
Mach Number of 2.0 With External Burning 



a 



• • 

gain in p^/p^ for m = 0.15 lb/sec is seen to be fairly in- 
sensitive to heat zone size and placement. 

Another set of similar calculations were made based on 
Q = 0.8 at w/r = 1.0. This value of Q corresponded to the 
design fuel flow rate of 0.6 pounds per second at 23,000 
feet altitude. The first three values of w/r, 1.0, 2.0 or 
3.0 would not compute because the external flow was driven 
to a sonic condition by the base pressure rise and the heat 
addition. This condition is analogous of thermal choking; 
however, in two dimensional flow both streamtube area change 
and heat addition drive the flow toward sonic conditions. 

The last three values of w/r = 4.0, 5.0 and 6.0, did not 
experience subsonic conditions in the external flow and 
their values of base pressure ratio are shown as the upper 
three curves in Fig. 22. The high value of p^/p^ = 1.91 
for the w/r = 4.0 curve was not expected. This represents 
a factor of six increase in base pressure over the two 
bottom curves in Fig. 22 with a factor of four increase in 
fuel flow. The combustion for this flight condition, Fig. 

23, was assumed to start at (z/r) = -1.2 and the initial 
compression was on the shear layer, downstream of the base, 
with full compression at z/r = +2.0. There is some error 
in the calculation introduced by neglecting vort icily and 
by the compression fan at the base, rather than the oblique 
shock. See Shaprio [Ref. 31] for a discussion of compres- 
sion fans. Assuming as much as a 1° error in the flow angle, 
far exceeding the estimate given in Section IV, the value of 
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P^/Pot remains at 1.8 and gives specific impulse of 230 
seconds. With w/r = 4.0 and z/r = 2.0, part of the com- 
pression due to combustion was applied along the shear 
layer downstream of the critical point. With z/r = 1.0 
more of the compression was applied upstream of the critical 
point, but the base pressure was less. The variation of 
P 5 /P 00 z / r indicates that the location of dp/ds (due 

to the external compression waves) along the shear layer 
may be an important consideration in achieving high values 
of base pressure with external burning. This hypothesis 
and the condition where the external flow with heat addi- 
tion becomes subsonic require additional study. 

Figure 23 provides some details regarding the optimum 
flow condition shown in Fig. 22 as point A. The heat zone 
and shear layer have been plotted in the lower half of 
Fig. 23. The boundary layer thickness has been included 
as well as the oblique shock wave which forms when p^/p^ 
is greater than one. For comparison, the shear layer 
without heat addition has been shown as a dashed curve in 
Fig. 23. Pressure ratios along the shear layer have been 
plotted in the upper half of Fig. 23; again, the dashed 
curve is without heat addition. In both cases, with and 
without heat, the shear layer pressure ratio shows a gradual 
increase with distance downstream from the base. The pres- 
sure ratio along the inner streamline of the heat zone 
(points a' through e') has also been plotted along with the 
shear layer pressure ratios (points a through f) . This 
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streamline experiences an abrupt increase in pressure ratio 
where it intersects the oblique shock wave at c ' ; then, its 
pressure ratio starts to decrease with a corresponding in- 
crease in Mach number, indicating that the flow along the 
inner streamline accelerates downstream of the oblique 
shock wave. The calculation along the shear layer was 
terminated before reaching a maximum pressure ratio due to 
the saddle point singularity, Fig. 4, since the straddle 
method employed in the computer program had established the 
required base pressure ratio within ± %% . 

Additional calculations with the computer program given 
in Appendix J would improve the accuracy of p^/p^ and would 
take values for the shear layer slightly closer to the 
critical point, but, since the critical point is a saddle 
point singularity, would not be expected to calculate through 
the critical point. Calculation "through" the critical 
point, using the straddle method, would require changes in 
the program to restrict the values allowed for dF/dw c and 
force the calculation across the singularity. Additional 
study is required to include this change and provide the 
capability of calculating the entire base flow field. 

To verify the existence of an oblique shock wave at 
the base, which was calculated in Fig. 23, one may consider 
the field of compression waves impinging on the shear layer 
as being concentrated into a single oblique shock wave im- 
pinging on the shear layer. This simplified condition, 
shown in Fig. 24, has been studied experimentally by Korst 
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Figure 24. Schematic Diagram of Resulting Base Flow When an External 
Shock Wave Impinges on the Shear Layer 



[Ref.* 36] and was found to create an oblique shock wave at 
the base. Thus, the calculated shear layer profile and 
pressure ratios in Fig. 23 arc in agreement with experimen- 
tal results from the simplified condition. 
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VI. CONCLUSIONS 



With the method of characteristics and the Crocco-Lees 
base flow model, one can accurately predict the base pres- 
sure ratio of a two dimensional or axisymmetric body with 
moderate external burning for various Mach and Reynolds 
numbers. The correct values for the two transition 
Reynolds numbers are essential for accurate base pressure 
ratio calculation of a particular body and should be ob- 
tained from experiment. The effect of projectile spin can 
be included approximately by dividing the transition 
Reynolds numbers with one plus the spin number. 

Net thrust can be produced by external burning. Plac- 
ing the combustion zone well forward of the base to give 
full compression on the body was found by experiment to 
produce a decrease in thrust. For maximum projectile 
thrust with supersonic external flow the computed heat zone 
placement was found to start 1.2 base radii upstream of the 
base and extend 4.0 base radii downstream. The zone started 
with an inner and outer streamline of 2.0 and 3.0 base radii 
respectively and received a fuel flow of 0.6 pounds per 
second of fuel rich propellant exhaust, ARC-168. It ap- 
peared higher values of base pressure would be possible with 
the same fuel flow rate by decreasing the length of the com- 
bustion zone still further. The flow field could not be 
computed with this heat distribution because of the resulting 
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subsonic flow and requires additional study to determine any 
further advantage or limits under these conditions. Experi- 
ments with actual combustion should be conducted to deter- 
mine the feasibility of the desired combustion placement 
and to measure the resulting base pressure. 

The rate of pressure increase along the shear layer due 
to the heat zone placement appeared to influence the com- 
puted value of base pressure ratio. This effect requires 
additional study to determine its significance in achieving 
higher values of base pressure ratio with external burning. 
Experimental studies are needed to determine any variation 
of the empirical F - k relation when external pressure waves 
impinge on the shear layer. 
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APPENDIX A: DERIVATION OP THE CHARACTERISTIC EQUATIONS WITH 

HEAT ADDITION 

Equations (1) through (4) govern steady compressible flow 
with heat addition [Ref. 18] . These equations neglect heat 
conduction, viscosity and body forces. In natural coordi- 
nates, Pig. 1, U = ws . The continuity equation, Eqn. (1), 
may be written 



where a = 0 for two dimensional planar flow and a = 1 for 
axisymmetric flow. The momentum equation, Eqn. (2), may be 
written for the two unit vector directions, s and n. 



Oi 

(2ur) pwAn = constant 



(1A) 




(2A) 



(3A) 



The energy equation, Eqn. (3), becomes 



3 s, 1 2-i 

K _ (h + _ „2) = q 




(4A) 



where the total enthalpy, h^ = h + ^w 2 , has been used. The 
thermodynamic relation for entropy, S, may be Avrritten 



TVS = Vh - - Vp 
P v 



(5A) 



t 3S 

3s 




(6A) 



t 3S 

3n 




(7A) 
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Using Eqn. (3A) into (7A) 



1 9w _ 99. = _ T 3S + _1 8h 0 
w 3n 3s w 2 " 3n w 2 3n 



(8A) 



Forming the logarithmic derivative of Eqn. (1A) in the s di- 
rection 



a 

r 



s i n 0 



+ 



1 9P + 
p 3 s 



where the relations 3r/3s 
have been used. Forming 
(4) in the s direction 



1 3w 39 n 

w 3s 3n U 



(9A) 



= sine and [1/An] [ 3An/3s] = 36/3n 
the logarithmic derivative of Eqn. 



I M = Hp + i II 

p3s p 3s T 3s 



Substituting Eqns . (2A) and (10A) into (9A) 



q • A 

— sm0 
r 



, — — 

1 w ' RT J 3s 



1 3T + 36 
T 3s 3n 



0 . 



( 10A) 



(HAD 



From Eqn. (4A) 



£L = 



w 




3 w 
3 



(12A) 



where the enthalpy has been written in terms of a constant 
Cp , h = c pT. Substituting Eqn. ( 1 2 A) into (11A) 



f 1 _ _w w n 3w_ 3_0 q asinO 

^ w RT c T J 3s 3n wc T r 

P P 



0. (15A) 



Combining terms, and using the Mach angle relation cot 2 y = 
M 2 - 1, Eqn. (13A) becomes 



cot 2 yi 3w 30_ asinQ _ q 

w 3s 3n r wc T 

P 



( 14 A) 
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Th6 Prandtl -Meyer Angle, v, is defined by the relation 



dv = dw . 

W 



Using Eqn. (15A) into ( 1 4 A) and multiplying by tan y 



9v . 90 . Jasin0 q f 

- tany = tany 4 - V > 

9s 9n ^ r wclj 



Similarly, Eqn. (15A) into (8A) becomes 
9v 90 T 9S 



1 9h 0 



tany 9n ‘9s w 2 9n + w 2 9n ’ 



Adding Eqns . (16A) and (17A) 



<-3^.* 9 w m , Iasin0 q 

C 3? tany to )(v ' e) “ tanu /- r 



wc T 
P 




Subtracting Eqn. (17A) from (16A) 




Now, in the characteristic coordinate system 

9 ( ) , 9 9 w . 

secy to = ( 5T + tmv to )( 5 

and 

9 ( ) ( 9 «. 9 w . 

secy - ( as - tany ^ )( ) . 



(ISA) 



( 16 A) 



(17A) 



(ISA) 



(19A) 



(20A) 



(21A) 
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With this, Eqns . (18A) and (19A) become 



(22A) 

3p - 0) = sinu f * sine g_ / _ cosy / 9S _ ^0 V. 

3ri ^ ^ r wc^T J w 2 ^ 1 3n 3n J 

and 

(2 3A) 

ai > ; °> - sinp ini . ♦ S°|E f r|i - 

3£ > r wc p / w 2 ( 3n 3n^ 

Eqns. (22A) and (23A) are the characteristic equations, with 
heat addition, governing the two dimensional planar, u = 0, 
or axisymmetric , a = 1, flow. These equations are dimensional, 
having the dimension of reciprocal length. 



78 



APPENDIX B: DERIVATION OF CONDITION FOR I RROTAT TONALITY WITH 

HEAT ADDITION 



The equations required for this derivation are momentum 
with no viscosity or body force, 



DU 

Dt 




(IB) 



where U is velocity, p is density, and p is static pressure; 
energy with heat addition, but without viscosity, body force, 
or conduction, 



D(h + P 2 U-U )_ n + 1 3p 
Dt 4 p 9 t 



(2B) 



where h is enthalpy, and q is time rate of heat addition per 
unit mass; the second law of thermodynamics 

TVS = Vh - ^ Vp ( 3B) 

where S is specific entropy; and the vector identity 

U-VU = V(%U-U) - UX(VXU). (4B) 

Combining Eqns. (IB), (3B) and (4B) 

— ■ + V(j U-U) - UX(VXU) = TVS - Vh . (5B) 

Introducing the definition for total enthalpy, h^ = h + !sU*U, 
into Eqn. (5B) , one obtains Croccos's Theorem 

TVS + UX(VXU) = Vh Q + . (6B) 

Let VXU = Si and take tlie curl of Eqn. (6B) , 

^ = vtxvs + fi-vu - nv-u (7B) 
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• ___ 

since V*ft = V*(VXU) = 0 and the curl of a gradient is zero. 
Eqn . (7R) describes the total rate of change of vorticity, 

ft . 

Now assume that ft = 0 prior to and during heat addition. 
Eqn. (7B) becomes 

= VTXVS = 0. (8B) 

Then, Eqn. (8B) is the condition which must be satisfied dur- 
ing heat addition if the flow is to remain irrotational ; 
that is, the temperature gradient must be parallel to the 
entropy gradient. Referring to Eqn. (3B) and assuming 

h = c T, where c is constant, the condition for irrotation- 
P P 

ality may be stated in terms of pressure; that is, the tem- 
perature gradient must he parallel to the pressure gradient. 

In the case of steady flow, the irrotationality condi- 
tion may be given in terms of heat per unit mass which has 
been added, Q. In Eqn. (2B) let q = U*VQ, and, using the 
total enthalpy h^ = h + ^U*U, 

U-V(h 0 - Q) = 0. (9B) 

Then along any streamline 

1>„ - Q - h (10B) 

where it is assumed that the upstream conditions, before heat 
is added, i.e. Q = 0, are uniform [Ref. 19]. With these as- 
sumptions one may take the gradient of Eqn. (10B) 

V(h Q - Q) = 0 (lib) 
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is constant. Using the steady form of Eqn. (OB), 



since h^ 



with ft = 0 , and Eqn. (11B) 



VS = - Vh Q = ~ VQ. 



(12B) 



Also, from Eqn. (11B) , one may write 



VQ = V (h + J^U-U) = c p VT + V(JsU«U). 



(13B) 



Combining Eqns . (12B) and (13B) in (8B) 




(14B) 



Since (VQ)X(VQ) = 0, Eqn. (14B) reduces to 



(VQ) XV (%U *U) = 0. 



(15B) 



Eqn. (15B) is the condition which must be satisfied to pre- 
serve irrotationality in steady flow with heat addition; that 
is, the gradient of heat added must be parallel to the gra- 
dient of kinetic energy. 

If the gradient of heat added is a small perturbation, as 
is the case of linearized flow, then the irrotationality con- 
dition in Eqn. (15B) is automatically satisfied. To show this 



is of second and third order perturbation. 

This condition may be applied to Eqns. (24A) and (25A) , 
Appendix A, and neglect the vorticity term 



let Q = h ro Q ' , where | Q ' | < < 1 , and U = U OT (i + g- ) , where U'*U' 



<<1, then, substituting in Eqn. (15B) , 



VQX V (%U • U) = h^VQ ' XV (U^i • U ' + ^U-U') 



(1GB) 
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thus- simplifying the method of characteristics calculation. 
First of all, if the heat zone is started in a region of 
uniform flow along a line normal to the streamlines, and if 
the heat is added uniformly in a zone bounded by two stream- 
lines, then the condition of Eqn. (15B) will be approximately 
satisfied within the zone since the gradient of heat added 
will be approximately parallel to the gradient of kinetic 
energy. Second, if the first condition is not satisfied 
(i.e., along the streamline boundaries of the heat zone where 
there is a jump in Q normal to the streamline), the vorticity 
may still be neglected if the heat added is a small perturba- 
tion, that is 



q = U • VQ = h^U • VQ ' . 



( 1 7 B ) 



Then the heat term used in Eqns . (24A) and (25A) must be 



wc T 
P 



h^wOQ'/as) 

wc T 
P 



• 9Q 1 

9s 



<< 1/r, 



(18B) 



This condition is satisfied for the heat zone boundary, where 
r^q/wCpT is taken as 0.171 to get - 4^° deflection of the 
inner heat zone streamline. 
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APPENDIX C: BRIEF DEVELOPMENT OF CROCCO-LEES BASE FLOW FORMU- 

LATION [Refs. 1,8,13] 

The mass flux, in, in the base flow is given by the inte- 
gral 

m = , P (2Trr) a pudr (1C) 

0 ^ 

where <5 is the viscous flow thickness, r is radius, p is den- 
sity, u is axial component of velocity, a equals zero for two 
dimensional planar flow and a equals one for axisymmetric 
flow. The momentum flux is 




( 2 C ) 



For the base flow, where there is no wall drag, the momentum 
equation becomes 



dl 

dx 



u 



dm 
e dx 



(tt6)“ 6 



, d J2 

dx 



( 3C) 



where u e is the external axial component of velocity, dp/dx 
is the static pressure gradient and, assuming small angles, 



dm #• a t 1 > n r dS * . 

33? = P e u e ( dx 6 o >• 



(4C) 



A major point of the Crocco-Lees theory is that 



- 0 

dx e 
where k is 

specified . 

mixing and 



= k ( 5 C ) 

a non dimensional mixing coefficient which must be 
The value for k is approximately 0.0053 for lamina 
0.03 for turbulent mixing. 
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* A mean velocity is defined by 

= I/m ( 6 C ) 

and a velocity ratio by 

k = u i/ u e ‘ ( 7C ) 

The mean density, p^, is defined by the equation 

m - ( Ti' 6 ) w Sp^u^ . ( 8 C ) 

It is assumed that the stagnation temperature is constant 
throughout the flow and that the lateral pressure gradient, 
dp/dr, is zero. Knowing the pressure, p = p , and p , a 
mean temperature, T^, can be defined 

T 1 = P e /p l R - ' (9C) 

Introducing the displacement thickness, 6*, and momentum 
thickness, 6**, of ordinary boundary layer theory 



<$* ( 1+ °0 = J 1 



(2r) a (l 



pu 

P u 
e e 



) dr = 6 



(1 + a) 



(IOC) 



m 



p U 7T 

K e e 



a 



and 



6 **(l+a) 





) dr 



= 6 



(1 + a) 



(1 + a) 



2 a 
p U IT 

K e e 



(11C) 



Noting that k = I/mu and using Eqns. (IOC) and (1JC), one 

G 

has 
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K = 



6 (l + ct) . 6 *( 1+a ) . 6 **( 1 + c 0 
6 (i + °0 . 6 *(i + a) 



(12C) 



On application of the Stewartson transformation [Refs. 8,13] 
one finds 6, 6*, and 6** in Eqn. (12C) are replaced by the 
incompressible quantities 6^, 6? and 6|*. Thus, within the 
restrictions of the Stewartson transformation, any compres- 
sible laminar base flow can be related to a corresponding 
incompressible laminar flow. It is assumed this also applie 
in the turbulent case. The ratio of mean temperature, T^, 
to stagnation temperature, T g , may be written 



r- f 

s 



y- 1 2 

2 1 



(13C) 



where 



K «( 1+a) 

l 

<$ (l + °0 . 6 *U + a) 
i i 



(14C) 



W 1 = U l /a s> anc ^ a s sta g nat 'i° n speed of sound. It has 

been found more convenient to work with the quantity F > 0, 



F + 1 = f/K 2 . (15C) 

Experimental data [Ref. 8] has shown a relationship between 
F and k which can be expressed by the polynomial 

K = 1.0 - 0.4F + 0 . 16F 2 - 0 . 0 11 8 5 F 3 . (1(3C) 



For F > 1.5, k = 0.72 and dK/dF = 0.0. 



From liqn . (8C) one may write an expression for mass flux 



m 



ma 



Y(Tr5) a 6 p e w e 

k[F + 1 - ^ w 2 ] 



(17C) 



Using this and Eqns . (3C) , (4C) and (5C) one may form a 
single ordinary differential equation 



dF 

dw 



rY + l 



w 



- 1) (1 - 



Y-l 

2 



w 2 ) + 

e J 




Eqn. (18C) has a saddle point singularity, [Refs. 1,8], also 
called the critical point, which, in integrated form, deter- 
mines the one unique solution describing the given base flow 
condition. In two dimensional planar flow without external 
heat addition, where 0 is known from w by the Prandtl- 
Meyer relation, one may start at the critical point, where 
both the numerator and denominator of Eqn. (ISC) are zero, 
and integrated toward the base [Ref. 8] . In flow with ex- 
ternal heat addition, and in axisymmetric flow, is not 
known as a function of w ; thus, in these cases, one must 
start integration of Eqn. ( 1 8 C ) at the base and concurrently 
solve the external flow field, by the method of character- 
istics or linearized theory, for the values of 0 along the 
shear layer. This process is one of trial and error described 
in section IV as the straddle method. 
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APPENDIX D: CALCULATION OF BOUNDARY LAYER MASS FLUX AT THE 

BASE 



The boundary layer is assumed to be adiabatic, have con- 
stant pressure across its thickness, 6, and have a velocity 

*j / 

profile (u/u ) = (y / 6 ) ' , where u is the external flow 

O G 

velocity and n is an integer [Ref. 24] ; 1 < n < 7, for 
laminar to turbulent flow. The boundary layer mass flux at 
the base is given by the integral 



m. 



- S 



x b* 6 



(2irr) a pudr 



(ID) 



where a is zero for 
axi symmetric flow, 
half-height. Using 
above assumptions, 




two dimensional planar 
p is density and r^ is 
the equation of state, 
Eqn. (ID) becomes 

r-r, 1/n 

(27rr) a u e (— ^ — — ) dr 

u 2 r-r, T”7 
T e ( b^2 /tl 

*T " 2c 1 6 J 

P 



flow, or one for 
the base radius or 
p = pRT, and the 



(2D) 



where p is the external pressure. Defining w as u /a , 

C 0 0 S 

where a s is the stagnation speed of sound, Eqn. (2D) becomes 



m b = m b a s = Y6p e w c (2u6) 



a 



r 

s 

r. 



r b + ^ 

r ^) 1/n (|) a C% (3D) 



(• 



i Y-l ..2 / _1 bx2/n 

1 — w c ( “6~ 3 
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Then 



m. 



M b = 



Y<Sp e w e (2TT6) 



p X 1/n (X+^-) a dx 

a J 1 - 



2 v 2 / n 



where x = (r - r,)/6 and B 2 = ^4- w 2 . 



2 e 

If a = 0 and x = z n , Eqn . (4D) becomes 

1 

‘dz 



2D 



■JV 



B 2 z 2 • 



0 



This may be integrated. For n odd, Eqn. (5D) gives 



(n - 1 ) / 2 



2D 



ODD 



= n f — hr ln ( J -r) - Yj — 

^ 2B n 1 1-B Z ^ [n+1 - 2 j ] B 



2 j 



(4D) 



(5D) 



where the summation is omitted if n = 1. For n even, Eqn. 
(5D) gives 



1 2D n 

ZU EVN 



2B 



If a = 1 and x = z 11 , Eqn. (4D) becomes 



n/2 

1 -j r 1+B , V 1 

n+ l (l - B 5 ' 4 i [„ * 1 - 2 _i 



(7D) 



2 i ] B 



2 j 



l AXI 



= n f T 



2n 



dz 



■B ? -z 2 



[> z n 

n j 1 - B 2 



dz 
B 2 z 2 



(8D) 



The second integral in Eqn. (8D) is the same as Eqn. (5D) , 
which is evaluated in Eqns. (6D) and (7D) . Integrating the 
first integral, Eqn. (8D) becomes 
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l axi 



^ 2B 2n + 



+ Ik i 

6 A 2D , 



n r l + B. 

T ln 1 -B^ 



V I t 

ftx [2(n-j)+l]B 2 V 



(9D) 



The results of Eqns . (6D) , (7D) and (9D) are used with Eqn . 

(4D) to evaluate m, or m, . This boundary layer mass 

2D d AXI 

flux at the base, before the corner expansion, m^ , is set 
equal to the mass flux after the corner expansion, m^ , ex- 
pressed in terms of the Crocco-Lees formulation, Appendix C. 



in, 



2D 



= Y«P w I 9n = m' 



yS ' PgWg 



/ e"e J '2D 11 b-,. .r-^, Y - 1 .2n 

2D k ' [ F 1 +1 - - - - - w ' ] 

i-t & 



(10D) 



Taking k' = 0.72 and the value of w^ assumed from the straddle 
method, one solves Eqn. (10D) for the value of F' which starts 
the integration along the shear layer. Similarly, in the 
axisymmetric case, 



m u = 2ry6 2 p w I. VT = m,* 

1 *e e AXI b 



7i y 6 ' 2 p ' w ' 

1 *e e 



AXI 



AXI 



k ' [ f ' + 1 - liiw; 2 ] 



(HD) 
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APPENDIX E: MASS FLUX CALCULATION ALONG THE SHEAR LAYER 



From Appendix C, Eqn. (17C) , the mass flux, m, may be 
written 



Y(iT6) a 6 p c w c 

m = = 

K [F+l-^~- w 2 ] 

1 2 e J 


(IE) 



where a = 0 for the two dimensional planar flow and a = 1 for 
axisymmetric flow. One may evaluate Eqn. (IE) at two adjacent 
points along the shear layer, m^ and and solve for the 

distance between these two points, dx , by using the additional 
relations [Ref. 1,8]. 



« = e * k 

dx 

and 


(2E) 


dm s /-v {s n cx 9 d 6 ^ 

37 " (2,t6:) V'e a s ( dx " 6 } • 


(3E) 



Assuming small angles, Eqn. (2E) states that the rate of 
change of the shear layer height or radius, d 6 /dx, is equal 
to the flow angle 0 plus some constant, k. Then the rate of 
change of mass flux, dm/dx, given by Eqn. ( 3 E) , is the mass 
entrained by the shear layer from the external flow. 

To compute m^ + ^, one increments F^ with AF^, say AF^ = 



0.1F. with the condition AF > 0.08, 
i - ’ 




F. , - = F. + AF. . 
x + 1 l l 

Then 


(4E) 


w = w + AF./(dF/dw ) 

i + 1 i 


(5E) 
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where dF/dw e is computed from Iiqn . (18C) . The value for K ^ + y 

is computed from an empirical relation using F^ + j, given by 
Eqn. (ICC), with = 0-72 if F^ + ^ > 1.5. The iscntropic 

relation between static and total pressure, in terms of y 
and w , is used to evaluate p^ at i and i+1. Average values 
between i and i+1 are used to evaluate d6/dx and dm/dx, where 
the isentropic relation between static and total density must 
be used with the latter. Since 

nr +1 = nr + (dm/dx) dx (4E) 

and 

6 i + 1 = 6. + (d6/dx) dx . (5E) 



One may solve directly for the value of dx , in the two dimen- 
sional planar case, which corresponds to the increment AF^ 



dx . 

l 



2D 



6 . 

i 




m i. + l , ( d6 dm 

6. J . 1 J / 1 dx ‘ dx J 
l + l 



(6E) 



For the axisymmetric case 



dx . 



1 AXI 



- 7 C dx. - T 1 44 dx? 



'2D 



6 . dx 

l 



L AXI 



(7E) 



Eqn. (7E) may be iterated several times, as long as dx. 



< 1, to determine the correct value of dx. 



x AXI 



l AXI 
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APPENDIX F: CALCULATION OF BOUNDARY. LAYER THICKNESS AND 

MIXING COEFFICIENT 



The laminar boundary layer thickness, 6j , was computed 
at point x from the compressible Blasius solution [Ref. 8] 

(IF) 

where x is the distance along the surface from the leading 



6 l = (5.0 + 1.2(Y-l)M 2 )x/(Re x )" 2 



edge and Re is the Reynolds number based on this distance. 

-A. 

If the base Reynolds number, Re^ , is less than the transi- 
tion Reynolds number, Re , then the boundary layer at the 
base is completely laminar and its thickness was computed 
from Eqn . (IF), with x taken as the body length. 

If the base Reynolds number, Re^ , is greater than the 
transition Reynolds number, Re^. , the boundary layer will 
be partly or completely turbulent. The turbulent thickness, 
6.p, was computed from [Ref. 8] 



<$ T 



^^(CfM/CfiHW**) 



( R e x ) 



(2F) 



Cpp is von Karman’s incompressible turbulent skin friction 
coefficient given by [Ref. 24] 



(0.242)/ (C £i ) 2 = 1 °gio ( ^ Rc x C fi ) 



(3F) 



and C^ is van Driest ' s compressible turbulent skin friction 
coefficient given by [Ref. 24] 



( 4 F ) 



0 



.242 (1-B Z )' 5 . -l n , , n „ , . i or, m „2, 

^ r 2 — sin B = lo glO^ Re x C fNp 1.261og 1Q (l-B ) 
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where B^ is (y-l)w^/2. Equations (3F) and (4F) must be 
solved by iteration. The quantity 6/6** (Appendix C) may 
be evaluated for the two dimensional l/7th profile using 
the same assumptions employed in Appendix D. Briefly 
+ 6 



6 ** = 



(P”) r i . ^dr 
(p u ) 1 u Jar 

r^ v e o J e 



(5F) 



ib +6 



V ■ c i ~t Lw ;p 



r- r 



(■ 



b } l/7 



r-r 



(- 



b } 2/7 



(6F) 






(dr/6) 



6 ** 

~6~ 



7 (1- B 2 ) (At In (- 



2B 



tt) 



1-B‘ 



V 1 r 

jti (8-2j)B^ 



- ~\ ln ( 

2B J 



1 + B- 
1 -B' 




1 

(9 - 2 j ) B 2 ^ 



). 



(7F) 



Since the actual transition process is more or less probabil- 
istic in nature, it was assumed that the transition from 
laminar to turbulent flow can be described by the integrated 
normal distribution. It was necessary to assume a turbulent 
transition Reynolds number, Re^ rt , for this purpose, where 
is approximately four times Re . The transition 
computation was made with the error function, erf, which is 
exactly related to the integrated normal distribution. The 
argument of the error function was formed in terms of base 
distance from the initial transition point, and was weighted 
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A numeri- 



so its value was -0.96 at Re , and 0.96 at Re 

* U 1 

cal span, S, of 4.0 takes the integrated normal distribution 
from 0.02 to 0.98. Formulating these conditions in terms of 
body radius, r, 

Z tr = (Rc tr /Re b ' i-OHL/r) (8F) 

is the transition point on the body measured from the base, 
and 



Z trt = ( Re trt/ Re b ' l-Om/r) (9F) 

is the turbulent transition point. Note that Z. . will lie 

1 trt 

downstream of the base if Re^ ^ > Re^ . Forming the average, 

Z trav " « Z tr * Z trt> ( 10F > 

is the distance from the base to the center of the transi- 
tion zone. If Z. is zero, the integral of the cumulative 
distribution must be 0.5. The value of erf at this point 
will be 0.0, and the cumulative distribution integral, 1^, 
is 



CD 



; ( 1 + erf) . 



The values of erf at the base, ERF^, range from 
+ 1.0. In terms of Z 

S Z 

ERF b = erf (— 



( 1 IF) 
1.0 to 



t rav 



2 2 (Z. . - ) 

v trt tr J 

Then, when Re^ > Re^ r , ^asE is given by 



(12F) 



'"'BASE = + ER V + - ERF b> 



(13F) 
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In ‘Eqn. (13F), ERF^ will be positive if Re^ is less than 

%(Rc,_ + Rc^ . ) , thus as Re, ranges from Re. to Re. . , the 

tr trt 7 ’ b b tr trt 

contribution of 6j will range from 0.99 to 0.01 and the con- 
tribution of 6^, will range from 0.01 to 0.99. The value for 
n used to compute the velocity profile at the base (see Ap- 
pendix D) is formed exactly as Eqn . ( 1 3 F ) 

"BASE = *(1 + ER V + ! U - ERI V U4F) 

since n ranges from 1 to 7 . 

When the transition zone lies along the shear layer, 
one must compute a transition value for the mixing coeffi- 
cient, k, at the point being considered, Xgj . The distance 



"BL 



is measured from the base in terms of r. To carry out 



this calculation, the value of erf at x^j , ERF PT , is written, 



BL 



as in Eqn. (12F) , in terms of Z 



ERF 



BL 



= crf(— 

2 2 



^trav X BL ^ 
< Z trt - Z tr) 



(15F) 



Note that Eqn. (15F) reduces to (12F) since x^j is zero at 
the base. Using Eqn. (15F) , the value for the coefficient, 
kpj , is formed in a similar manner from the laminar and 
turbulent mixing coefficients, kj and k,j, , 

k BL = ^ k L (1 + ERF BL ) + ^ k T ^ 1 ‘ E r F bl ) . (16F) 

In Eqn. (16F) , the value of ERF ni will be positive if Z 

15 L) 1 1 n v 

is downstream of x^j , and vice versa. 
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APPENDIX G: DERIVATION OE LINEARIZED STEADY PLOW EQUATION 

WITH HEAT ADDITION 



The equations, (1) - (4), are required for this deriva- 



tion . 

V- (pU) = 0 (1G) 

pU-ViT = - Vp (2G) 

U • V (h + %iJ‘U) = q ( 3G) 

p = pRT (4G) 

Forming the directional derivative, U*V, of Eqn . (4G) 

U-Vp = pRU-VT + RTU-Vp. ( 5G) 



Substituting Eqn, (1G) into (5G) 
pRU • VT = U-Vp + pRTV • U 

Substituting U' of Eqn. (2G) in (6G) 
RU • VT = - U-V(%U-U) + RTV • 0" 



(6G) 



(7G) 



Assuming is constant, using h = c^T and c^ = yR/(y-l), 
substitute y/(y-l) times Eqn. (7G) into (3G) 



V-U - U-V^U-U) = q/ c T. 

a Z p 

Let T = T m ( 1 + T'), q = h m q ' , U = UjT + V4> 
With this, Eqn. (8G) becomes 
(U T + V<f>) • V 



C 8 G ) 



V"<|> - -o 

a; (1 + T') 



(bU* + + W-V<fr) 



(9G) 

q'/Cl + T') 
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Retaining only first order terms, Eqn. (9G) becomes 



V 2 <j> - M 2 <j> = q’ . 

T co Y xx M 

This is the linearized steady flow equation, with heat 
tion for irrotation subsonic or supersonic flow. 



(10G) 
addi - 
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APPENDIX H: DETERMINATION OF THE MACH ANGLE 



This appendix derives a numerical procedure for deter- 
mining the Mach angle, y, for a given value of the Prandtl- 
Meyer angle, v 

The Prandtl -Meyer angle may be written [Ref. 18] 

v = — tan ^(g(M 2 -l)^) - tan ^(M 2 -l)^ (1H) 

g 

where M is the Mach Number, 1 < M < 00 , 

g = C(Y-l)/ (Y + l) ) ?S (2H) 

and y is the specific heat ratio. Let 

y + e = ir/2 (3H) 

where y is the Mach angle defined by 

sin y = 1/M. (4H) 

From Eqns . (3H) and (4H) 

(M 2 -l) % = tane. (5H) 

Using equation (5H) , equation (1H) may be written 

tang(v+e) = gtane. (6H) 

In Eqn . (6H) , the relation between v and e is still implicit. 

For y = 5/3, g = h and Eqn. (6H) can be reduced to the ex- 
plicit form 

tan 3 y = tan^- . (7H) 

Thus, using Eqns. (311) and (711), the simple relation may be 
written when y = 5/3 
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y = 2. - 2 tan 1 ( (tan^-) 1//3 ) . (8H) 

This analytical result [Ref. 37] was first given in terms 
of M, rather than y, along with a more complicated expres- 
sion when y = 5/4. There are no analytical expressions, 

V* = y(v,y), for other values of y, consequently, y must be 
determined from tables, or numerically, in all other cases. 

Equation (6H) may be written in a form for iteration 
which always converges to the correct value of e 

e i+l = tan_1 (| tan g( v+e i) ) (9H) 

for 0 < e. < tt/2 and 0 < v < v . The value v is given 
x - - max max & 

by Eqn. (1H) for M infinite 

v max “ ? ((CY + 1)/(Y-1))' 5 - 1). (10H) 

One may start the iteration of Eqn. (9H) with an approxima- 
tion for given by Eqn. (7H) 

= 2tan ^ ( ( tan^-) ^^) . (lid) 

This approximation fails for large values of v since it may 
give a value > tt/2. When this occurs, set = tt/2 to 
start the iteration with Eqn. (9H) . 

While Eqn. (911) will always converge to the correct value 
of e, its convergence is slow and, after two or three itera- 
tion, should be replaced by a method utilizing the derivatives. 
In effect, Eqn. (911) is used to find an approximation of e 
which lies in the range of convergence of a faster method. 
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One method which provides rapid convergence may be written 
from Eqn. (6H) by using a Taylor's series expansion for the 
separate sides about and solving for Ac. Let 



f 1 = tang(v+e i ) 
i 

f 2 = gtane^ 
i 

then 



Ac 



j + l 





(f{” - fj") + ---^f^ - f{ ) 

i i i i 



( 1211 ) 

(13H) 



(14H) 



where Ae^ = 0 starts the iteration in Eqn. (14H) . The advan- 
tage of this form lies in the fact that the derivatives may 

be written in terms of the values of f. and f_ , i.e. fi 

1 Z i l 

and f£ involve the sec 2 x, which can be written as 1 + tan 2 x, 

l 

etc. Thus the derivatives may be included with only two 
numerical evaluations by the computer of the tangent functions 
in Eqns . (12H) and (13H) . Depending on the number of deri- 

vatives included, two or three iterations with Eqn. (14H) 
should be sufficient. Then, of course. 



e. Al = e. + Ae.., (15H) 

i+lij+1 v ' 

is used to evaluate f. and , etc. For y = 1.4 and 

i+l l + l 

including the fourth derivative, the above procedure will 



provide twelve place accuracy for e in two iterations. With 
the required value for c, y can be written using Eqn. (311). 
This method was used in the computer program discussed in 



Appendix J. 
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APPENDIX I: SPINNING CYLINDER DESIGN 



The spinning cylinder was designed to be compatible with 
the coaxial free jet fixtures, nozzles and spacers used for 
the nonspinning investigation [Ref. 29]. The spinning cyl- 
inder and its support assembly had the same external dimen- 
sions as the stationary cylinder design, 1.282-inch diameter 
and 9.95-inch length, to insure similar test results without 
spin. The cylinder was driven by an internal impulse tur- 
bine. The design speed of 48,000 revolutions per minute, to 
match the spin number of 0.1595, was not achieved. This was 
due in part to physical interference during the fabrication 
process which reduced the tangential component of the turbine 
jet to less than one-half the desired value. The maximum 
speed achieved, without nozzle flow, was 18,000 revolutions 
per minute. 

The spinning cylinder was made of Plexiglas, to reduce 
rotating mass, and was glued to a metal endcap which formed 
the cylinder base. The endcap was driven through a shaft 
which screwed into the internal turbine. The shaft also 
carried thrust loads, from an internal air bearing surface 
on the turbine face, to react with pressure and friction 
forces due to the nozzle flow. This part of the design was 
deficient in that the loads due to nozzle flow caused the 
turbine face to rub and stopped the cylinder. Time did not 
permit modification to correct this deficiency. 
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The spinning cylinder was positioned radially by 24 air 
bearing pads provided along the outer surface of the support 
cylinder. High pressure air was supplied to these pads 
through orifices formed by drilling into supply tubes sol- 
dered into the walls of the support cylinder, Fig. 11. The 
bearing air exhausted through ports to the inside of the 
support cylinder where it mixed with turbine exhaust and 
flowed out the bottom. 

The base pressure on the spinning cylinder was measured 
by a total pressure probe, positioned on centerline next to 
the base, and a mercury manometer. The cylinder speed was 
determined through a photo diode and a pulse counter. The 
photo diode was illuminated by a laser beam reflecting from 
the metal endcap and half of the endcap was painted black to 
generate one pulse per revolution. 



Figure II. Disassembled View of Spinning Cylinder. 
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APPENDIX J: COMPUTER PROGRAM 



The FORTRAN computer program, included at the end of 
this appendix, required 150K for compile, 110K for link and 
80K for execution. It was written to compute base pressure 
ratio, within approximately 0.5 percent, by applying the 
equations and methods discussed in this thesis. In view of 
the empirical information and approximations employed in 
many of the equations, computation to greater accuracy did 
not appear warranted. Using an IBM 360 computer, computation 
of the base pressure ratio for one flight condition required 
approximately two minutes for the fully turbulent shear 
layer without heat addition. The addition of heat upstream 
of the base and low Reynolds number, which involved a long 
laminar shear layer, increased the size of the characteristic 
mesh and resulted in computation times as high as ten minutes. 

The flight condition was specified before making a com- 
puter run. This included choosing the Mach number, card 80, 
and the base Reynolds number, card 140. The two no spin 
transition Reynolds numbers, cards 120 and 130, were choosen 
for the particular body to provide agreement with available 
experimental data. The ratio of body radius to body length 
was specified on card 150. The spin number, rw/U, was given 
on card 131. The ratio of specific heats for air was set at 
1.4 on card 200. Card 180 and 190 were interchanged for two 
dimensional planar flow. For flow without heat addition, 
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car‘d 90 was set to zero, card 60 to 1STRT = 90, and card 100 
to ZSTRT = 20.0. 

Several additional considerations were required to com- 
pute the flight condition with external burning. The start- 
ing values for the radii of the heat zone inner and outer 
streamlines were specified on cards 420 and 430. The axial 
positions to start and stop heat addition were given on 
cards 100 and 110, with reference to the base. The value 
for Q, card 90, was specified to obtain a desired - 4h° de- 
flection of the inner heat zone streamline, however, this 
value could have been specified by assuming various fuel 
flovA rates and heating values. Finally, the value for ISTRT, 
card 60, was set low enough (minimum value of 1) to place 
the Mach line boundary condition well upstream of the outer 
compression waves generated by the heat zone. The values of 
selected quantities were printed, to monitor the progress of 
the computation, in addition to the required values for the 
base pressure ratio. 

A simplified flow diagram is shown in Fig. 1 J . It is 
intended to show the basic sections of the computer program 
and computational details have been omitted. 

In the computer program the group of IF STATEMENTS which 
examine the mesh point for heat addition can introduce a 
slight error at two corners of the heat zone. The amount of 
heat assigned between two points is proportional to the ele- 
ment of length along the characteristic in the heat zone. 

The IF condition is checked in radial position against the 
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Figure 1J. Simplified Flow Diagram for External Burning Assisted 
Projectile (FBAP) Computer Program. 
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Inner and outer streamlines of the heat zone, and in axial 
position against the start and end of the heat zone. Thus 
a right running element at the lower left or upper right 
corner of the heat zone could be outside the zone, but 
still pass the IF condition tests and be assigned % value 
for heating. This error is small with a small characteristic 
mesh and is neglected. The same applies to a left running 
element at the upper left or lower right corner of the heat 
zone . 
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